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sider a cosmological model where the cosmological constant Ao is replaced by a cosmological 
term A(t) which is allowed to vary in time. More specifically, we are considering that this dark 
energy term interacts with dark matter through the phenomenological decay law p\ = —Qp\. 
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(SNela, GRB, CMB, BAO, OHD), emphasizing the case where n = 3/2. This case is the 
only one where the late-time value for the ratio of dark energy density and matter energy 
density p\/p m is constant, which could provide an interesting explanation to the coincidence 
problem. We obtain strong limits on the model parameters which however exclude the region 
where the coincidence or the cosmological constant problems are significantly ameliorated. 
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1 Introduction 

It is now widely accepted by the cosmological community that our Universe is currently 
experiencing a phase of accelerated expansion. The first significant evidences came from 
the distance measurements of type la supernovae (SNela) [1-3] (see [4, 5] for an update) 
which has been later confirmed, among others, by the measurement of the anisotropies of 
the cosmic microwave background (CMB) spectrum by the Wilkinson Microwave Anisotropy 
Probe (WMAP) [6-9] and by the measurement of the baryon acoustic oscillations (BAO) in 
the Sloan Sky Digital Survey (SDSS) luminous galaxy sample [10, 11]. 

Normal matter satisfies the strong energy condition (w = p/p > 0) and cannot drive the 
accelerated expansion of the Universe. Hence, the recourse to a cosmic fluid, known as dark 
energy, characterized by a negative pressure (w x < 0), is usually made to explain the current 
accelerating expansion of the Universe. Alternatively to dark energy, other approaches such 
as modified gravity [12-14] and inhomogeneous cosmology [15, 16], have been proposed to 
explain the (apparent) cosmic acceleration. 

The simplest form of dark energy is provided by a cosmological constant, which has 
led to the development of the ACDM model. In this model, the Universe is composed of, 
in addition to ordinary matter (radiation, w r = 1/3, baryon, Wj, = 0), a pressureless cold 
dark matter fluid {w c = 0) and a cosmological constant Ao (w x = = —1). The model 
provides a reasonably good fit to the current cosmological data but is however is plagued 
by two serious theoretical difficulties. One of them is the cosmological constant problem; the 
observed value of the cosmological constant is ~120 orders of magnitude smaller than what 
is expected from theoric computations [17]. The other one relies on the observation that 
the values of the matter energy density (p m = pb + p c ) and of the dark energy density are 
currently of the same order of magnitude. However, since the former is diluted proportionally 
to the volume of the Universe as it is expanding (p m oc a -3 ) while the later remains constant 
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(pa = PA )) the period of time during which p m /p\ ~ 0(1) is predicted to be very narrow 
in comparison to the whole Universe history. To currently lie in this narrow period, a fine 
tuning of the initial conditions of the model is needed. This is the essence of the coincidence 
problem. 

A possible approach to alleviate the cosmological constant problem consists in replacing 
the cosmological constant Ao by a cosmological term A(i) whose value varies over time. This 
can also help to decrease the severity of the coincidence problem if, for instance, the matter 
dilution caused by the cosmic expansion is compensated (at least partially) by the matter 
produced through the dark energy decay. A cosmological model with a dynamical A has 
been proposed as early as 1933 [18] and the idea regained attention in the 1980s [19-23] 
following the publication of studies [24-29] suggesting that the dynamical effects of quantum 
fields in de Sitter space-time would lead to a decay of the effective value of the cosmological 
constant (see refs. [30-32] for an update). In the absence of a complete theory of quantum 
vacuum in curved spacetime, several decay laws, mostly based on empirical arguments, have 
been proposed in the literature. A review of the first A(t)-models developed may be found 
in ref. [33] and some more recent examples in refs. [34-47]. 

In this study, we will consider the phenomenological decay law 

Pa = -QpI (1.1) 

where Q is a parameter which can be constrained from observations and n an index charac- 
terizing different models. This law has been proposed in ref. [19] to study the decay of dark 
energy into radiation, but it was shown that the constraints on this process are too tight to 
allow p\ to decay from an initial large value to that observed today. The possibility that 
dark energy decays into ordinary matter is ruled out because the annihilation of matter and 
anti-matter would produce a 7-ray background in excess of observed level [20]. In this paper, 
we will then only consider the interaction between dark energy and dark matter. Although, 
only the decay of dark energy into dark matter [Q > 0) is relevant to possibly explain the 
cosmological constant problem and the coincidence problem, we will also consider the decay 
of dark matter into dark energy {Q < 0) for the sake of completeness. 

The model with n = 1 has already been studied as a peculiar case of two more general 
models, one where the decay depends also on the cold dark matter density (p\ = QaPA+QcPc) 
[34, 35] and one where the equation of state parameter of dark energy is not fixed to — 1 
{p x = —Qpx) [36]. Another generalization, involving this time a generic value of n, has been 
proposed in ref. [37] {p\ = —Qp^p^ 1 ); however the case where eq. (1.1) is recovered (n m = 0) 
did not have been considered. Although all these generalizations would be interesting to 
incorporate in our study, the resulting analysis would be significantly complicated. We will 
thus restrict our attention in this paper to the simple case of eq. (1.1). 

The model with n = 3/2 is particularly interesting with regard to the coincidence 
problem since it is the only one where the ratio p m /pA is constant at late time. For this 
reason, we will give it with a particular attention in the remainder of this article. It is to be 
noticed that solutions found in refs. [38-42] where the ratio p m /pA is constant for all time 
and where w x = wa = — 1 actually correspond to a subclass of solutions of this model (see 
section 2.1 for more details). 

2 Model 

We consider a model where dark energy and dark matter interact with an energy transfer 
rate Qa = —Qc = ~QPa m a spatially flat Friedmann-Robertson- Walker (FRW) spacetime. 
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The sign of the constant Q sets the direction of the transfer. The interaction is a decay from 
dark energy to dark matter (A — > DM) for Q > and a decay from dark matter to dark 
energy (DM— > A) for Q < 0. The energy conservation equations for radiation (w r = 1/3), 
total matter (dark and baryonic, w m = 0) and dark energy (wa = — 1) read respectively 

p r = -4H Pr , (2.1) 

p m = ~3Hp m + Qpl, (2.2) 
PA = -Qpl, (2.3) 

and the Friedmann equation takes the usual form 

3H 2 



Pr + Pm+ PA- (2.4) 



If we set to = (the variables with the subscript refer to their present value), the 
solution to eq. (2.3) is given by 




PA(t) = < 
where t n is defined as 

-O-i) 

'"^•^-(^W < 2 - 6) 

In the limit where Q approaches 0, the time t n becomes infinite and we recover the ACDM 
model (pa = PA )- For the decay of dark energy (Q > 0), the time t n is situated in the future 
(t n > 0) if n < 1 and in the past (t n < 0) if n > 1. The converse is true for the decay of 
dark matter (Q < 0). Depending on its sign, t n constitutes, or a past bound, either a future 
bound on the validity of eq. (2.5). Indeed, for n > 1 the model clearly breaks down when 
t approaches t n since p\ becomes infinite. For n < 1, the dark energy density is p\ = at 
t = t n . Hence for Q > 0, t n corresponds to the latest time when the interaction could be 
physically meaningful since the dark energy would have completely decayed in dark matter 
at this point. Similarly, for Q < 0, t n corresponds to the earliest time that the interaction 
could be meaningful. 

As mentioned in the introduction, the case where n = 3/2 is the only one which leads to 

a constant ratio of p m /pA at late times and could therefore provides an elegant explanation 

to the coincidence problem. Indeed, for n ^ 1 and \t/t n \ S> 1, we see from eq. (2.5) that p\ 

i 

becomes proportional to t . Moreover, at late times p m is by hypothesis proportional to 

PA and the radiation can be neglected in the Friedmann equation, hence the Hubble parameter 

i 

H becomes proportional to t 2n ~ 2 . It is easy to verify from eq. (2.2) that the only value 
of n consistent with these expressions is 3/2 (otherwise, the term Hp m would not have the 
same time-dependence as the other two). Because of this remarkable feature, we will consider 
in details the model with n = 3/2 and find analytical solutions in order to have a deeper 
understanding of its characteristics. For that, we will divide the cosmological evolution into 
two eras: an early one dominated by radiation (RD-era), and a late one dominated by matter 
and dark energy (MAD-era). 

For all values of n (including n = 3/2), we will use numerical solutions to constrain the 
models parameters. In that case, it will be more convenient to express the energy conservation 
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Figure 1. The curve fl\ = shown in the 

plane il\ a — Q. The parameter Q is expressed here 
in term of the dimensionless quantity Q/G 1 * , where 
G is the Newton's constant. The physical solutions 
(k < 0) are situated below the curve f2A = ^A^ ■ 



Figure 2. The ratio PkJpk obtained from 
cq. (2.19) and presented for different values of ^Ao 
(indicated on the figure) as a function of the dif- 
ference between Qa^ and Oa - The dashed line 
represents the ratio between the theoretical value 
and the observed value of dark energy density ob- 
tained in the context of the usual ACDM model 
(~ 10 120 ). 



equations (2.1-2.3) as functions of the redshift (corresponding to data) 



dp r _ 4 Pr 
dz 



dp r . 



(! + *)■ 

Pm 



dz 
dpA 



Q- 



Pa 



-Q 



(i + z) ^H(i + Z y 

Pa 



dz " B H(l + z) ' 

2.1 Matter— Dark energy dominated era (MAD-era): n = 3/2 



(2.7) 
(2.8) 
(2.9) 



During the MAD-era, we assume that p r <C p m , pa- We can then write the Friedman equation 

as 

3H 2 



8vrG 



Pa- 



(2.10) 



In order to have a solution for p m , we must solve for H first. Using the preceding equation 
to replace p m and p m in eq. (2.2) leads to the equation 



AttG 8ttG 
The solution for n = 3/2 is given by 



Hit) = S Q 



(2.11) 





l 




PA 


2 








_ \8ttg) _ 







Mpa/pa ) 2 - 1 
k{pA/pA )~^ + I, 



(2.12) 
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(a) (b) (e) 




Figure 3. Examples of the evolution of the dark energy density parameter 17a (thick line), the 
matter density parameter fi m (thin line) and the modulus of the scale factor (|a| — for t < ti, a S I) 
in absence of radiation and baryons for an expanding Universe (Ho > 0) with (a) Oa > ^A^ , (b) 
fiAo < ^A=o < 1 an d (c) f^Aoo > 1- The dashed lines correspond to the solutions that we would obtain 
by letting p m become negative and the shaded areas to the mathematical solutions for times prior 
to the singularity. We have considered that the interaction starts at t on = t w for (a) and stops at 
toff — t w for (c). For (a), we can notice the absence of an initial singularity and of an early phase 
dominated by matter 



where A; is a dimensionless integration constant, Sq = sign Q and /3 = Sq^/I + 96irG/Q 2 . 
The scale factor a can be obtained by integrating H, 



a = aoexp 



/**) 



a 



=1+1 / „£\ 2 /3' 
p A b (l + k(p A /pA ) 2 



Sq 



In order to have a = 1 at t = to, the value of the constant clq is set to 



a 



Pa 6 (1 + *) 



2/3 



1-Sq 



(2.13) 



(2.14) 



The solutions obtained can be divided into two classes depending on the sign of the 
constant k. If we only consider dark energy and dark matter, it can be shown that the 
solutions with k > are characterized by the absence of an initial singularity and by the 
absence of an early phase dominated by matter (which is needed for structures formation). 
The converse is true for k < 0, i.e. the solutions are characterized by the presence of an 
initial singularity and by the presence of an early phase dominated by matter. 

To find an expression for k, the initial conditions [p\ , Q, Hq] could be used, however, it 
will be more convenient to replace the first two parameters by the current and the late-time 
value (supposed to be constant for n = 3/2) of the dark energy density parameter We 
can rearrange eq. (2.12) in order to have 



M ' Ptot Ml) 

\8ttG J 



1 



1 + 



k(pA/pA ) 



1 



k(pA/pA ) 2 + 1 



(2.15) 
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In the limit where t — > oo, p\ becomes as expected and Oa approaches a constant value, 
JIa —t- . The parameter f^A^ depends only on the value of Q and is defined as 



^A 



til 

/3 + 1' 



(2.16) 



We have to keep in mind that £l A Q truly represents the late time value of ^l\(t) only for the 
decay of dark energy. For the decay of dark matter, we have shown that the model breaks 
down in the future at t = t n , hence before reaching the "late time" (actually, as we will see, 
the model breaks down even before t n ). In this case, ^a^ has to be considered only as a 
parameter related to the strength of the interaction. Setting t = to in eq. (2.15) and solving 
for k yields 



ft! 



(2.17) 



where Sh = sign Hq. For an expanding Universe (Hq > 1), the sign of k depends only 
on whether ^a^, is greater or smaller than Qa - Since Q,\ and n A ° are D °th included in 
the interval [0, 1], A; is necessarily negative for the decay of dark matter but could be either 
positive (if Qa > ^A^) or negative (if Qa < ^A^) for the decay of dark energy. The 
solutions with k = correspond to a situation where the ratio p m /pA is constant for all time 
(at least as long as radiation is neglected) and which is a peculiar case of the solutions found 
in refs. [38-42]. The curve = ^a^ is shown in the plane f^Ao — Q in figure 1. For the 
decay of dark energy, the relevant range of values for Q is roughly situated between 10~ 2 G2 
and 10 2 G2. Below the lower bound, the interaction becomes insignificant and we recover 
approximately the ACDM model (f^Aoo ~ 1) ; whereas above the upper bound, Qa remains 
close to zero for all time. 

From eq. (2.13) it should be obvious that only the solutions with k < feature a 
singularity and that it occurs at the time 



1 



(2.18) 



For Hq > 0, this singularity occurs in the past (ij < to) and could be considered as the initial 
time of the Universe. Evaluating eq. (2.5) at t = tj, we find an expression for the ratio of the 
initial and the current value of the dark energy density 



^ = (-*)* 

PAo 



(2.19) 



The case where J7a ~ ^Aoo could in principle provide an explanation for both the 
coincidence and the cosmological constant problems. Indeed, when radiation is neglected, we 
have Pm/pA = 1/^A — 1) hence the condition £Ia ~ ^A^ implies that the current value of 
the ratio p m /pA is now typical of a large period of time and there is no more "coincidence" 
problem. Moreover, if Qa ~ ^A^; the ratio PAi/pA can be arbitrarily large. However, if the 
difference between the theoretical value of pAi and the value of pa obtained for the A(t)CDM 
model remains of several orders of magnitude, we see from figure 2 that a fine tuning of the 
quantity SIa^ — ^a would be needed to explain it. 
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To complete the analysis of the MAD-era, we must come back to the problem encoun- 
tered in the future for the decay of dark matter. It easy to verify that at t n , the density 
parameter of dark energy is greater than 1, which constitutes a violation of the null energy 
condition (w < —1). Setting Q\ = 1 in eq. (2.15) and solving for t, we actually find that the 
null energy condition starts to be violated at 



1 _ ( lp 



t n , (2.20) 



where 

i 

k w [n Ao ,n Aoo ,H ] = 5fc5 Q fl A o + "A 00 (2 21) 

SkSQ^i^ - 1 

and Sk is defined as the sign of the constant k. The density parameter fl\ could become 
greater than 1 only if we let p m become negative after t w . To prevent this problem, when 
the dark matter has completely decayed in dark energy (p c = 0), we should consider that 
there is no more interaction and the dark energy density becomes a constant. If we neglect 
the baryon density, the interaction stops at t g = t w and the late time(i > t ff) value of p\ 
is then given by 

PA» = (-£)%A - (2-22) 

The solutions for SI a, Q m = 1 — and \a\ are shown in figure 3 for VL\ > ^Aoo) 
Oa < ^Aoo < 1 and ^a^ > 1- As we can see, the null energy condition is also violated for 
the decay of dark energy, but this time the violation occurs in the past. The presence of the 
singularity prevents this problem (ij > t w ) for k < but not for k > 0. We could simply 
consider that p\ is constant before t on = t w . However, the reason why the dark energy 
should start to decay at t w would be unclear. It would also be possible to consider that 
decay starts after t w (when p m is non-zero), which would lead to a solution with an initial 
singularity and an early phase dominated by matter. That could be the case for instance if, 
for some reason, the dark energy decay starts only when the dark matter becomes sufficiently 
diluted. However, we do not have considered any threshold of this kind for the case k < 0. 
In order to be consistent and to keep our model as simple as possible, we will not introduce 
any threshold, but rather consider that the model is not valid for k > 0. 

2.2 Radiation dominated era (RD-era): n = 3/2 

Until now, the radiation has been neglected. This approximation clearly does not hold at 
early times. Solving eq. (2.1) leads to the usual result (which holds for all time) 

p r = p ro a~ 4 . (2.23) 

During the RD-era, we can neglect the dark energy and the matter densities (p r ^> p m ,p\). 
Using the Friedmann equation, we can easily find the Hubble parameter 

H = 2H ro a~ 2 . (2.24) 



where H ro = \J (87rG/3)p ro and integrate it to obtain an expression for the scale factor 

a = [2H ro (t - t eq ) + a 2 eq ]K (2.25) 
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Figure 4. Evolution of dark energy (Ha, thick line), matter (fi m , thine line) and radiation (f2 r , dashed 
line) density parameter for different values of SIa^ and with H,\ = 0.73 and Ha = 70 km s _1 Mpc _ . 
The values of f2& and Q ro are set using eq. (3.12) and eq. (3.13). The time has been redefined setting 
U = 0, hence the age of the Universe is given by to. 



The value of the integration constant has been set in order to have a eq = a,MAD(t eq ) = 
a RD(teq), where t eq is defined as the time where p r = p\ + p m . The initial singularity now 
occurs at 

U= teq (2.26) 

We cannot find an analytical expression for t eq ; instead, we have to make the approximation 
that the solutions found for a and H during the MAD-era hold up to t eq and solve numerically 
the equation 

M + P^MAD\teq)i 

1H 2 (t ) (2-27) 

f J ro a MAD\ L eq) ~ g^Q 

To find an expression for the matter density p m during the RD-era, we can use the expressions 
found for p\ and for H (eqs. (2.5) and (2.24)) to solve the matter conservation equation 
(eq. (2.2)). Setting the value of the integration constant in order to have p meq = Pm RD (t e q) = 
Pm M AD(teq), we get 

p m = p meq (j-J + [*!(*) + F 2 (t)} a~ 3 , (2.28) 

where F\ is given by 

*i = {p^Ac - PAa 3 ) + 3H ro \t n \pl o (pl eq a eq - pfa) (2.29) 
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and F2 by 



F = 



a 



arctan ( — 
a 



arctan 



/9A , with a = J\ 2H ro (t eq - t 



In the limit where Q approaches 0, F\ and F2 are reduced to F\ 



-F2 — 3H ro \t n \(a, 



eq 



(2.30) 



so the second term in eq. (2.28) vanishes and we recover the expected expression for the 
ACDM model. 

As a consequence of the approximations made for the MAD-era and for the RD-era, 
Ptot 7^ 3H 2 / (8ttG) in the vicinity of t eq . To evaluate ptot (and Oj) we should instead explicitly 
add pa, p m and p r together. The evolution of the density parameters are shown in figure 
4 for different values of . For smaller values of ^a^ , the age of the Universe is larger. 
There is also an inverse relationship relating f^Aoo and the duration of the RD-era and of 
the AD-era, whereas the duration of the MD-era decreases as f^Aco becomes smaller. The 
variation of these time intervals becomes increasingly important as approaches £l\ , 

leading to solutions significantly different from what we would have obtained in the ACDM 
model, even for t < to- We can then expect that these solutions will be disfavoured by the 
observational constraints. 



3 Observational Constraints 

In this section, we explain our methodology for using the currently available data to constrain 
the three free parameters (S1a ? ^Aoo , Hq) of the model. Following what was done in refs. [43- 
46] to constrain cosmological models with interacting dark sectors, we will consider the tests 
described below involving the distance modulus p of type la supernova (SNela) and gamma- 
ray bursts (GRB), the baryon acoustic oscillation (BAO), the cosmic microwave background 
(CMB) and the observational Hubble rate (OHD). In refs. [43-46], the authors also use the 
gas mass fractions in galaxy clusters as inferred from x-ray data to constrain their models, 
however since this method involves many uncertain parameters, we do not consider it here. 
The best fit follows from minimizing the sum of the x 2 of each data set 

Xtot = XI + XOHD + X 2 BAO + XCMB- C 3 - 1 ) 

3.1 Distance modulus p of SNela and GRB 

The distance modulus is the difference between the apparent magnitude m and the absolute 
magnitude M of an astronomical object. Its theoretical value for a flat Universe is given by 

p th (z) = 5 log 10 + 42.38, (3.2) 

where h = i?o/(100 km s _1 Mpc _1 ) and the Hubble free luminosity distance Dl is defined as 

D L = H (l + z) I —. (3.3) 

The best fit is obtained by minimizing the x 2 function 

x 2 [^o^A OC ^o]=E [ ^ (Zt) ; 2 ^ (Zt)]2 > m 

i 1 
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z 


[58] 


0.1 


0.17 


0.24 [59] 


0.27 


0.34 [59] 


0.4 


0.43 [59] 


H(z) (km s^Mpc" 1 ) 


74.2 


69 


83 


79.69 


77 


83.80 


95 


86.45 


l-<7 uncertainty 


±3.6 


±12 


±8 


±3.61 


±14 


±4.55 


±17 


±4.96 


z 


0.48 


0.88 


0.9 


1.3 


1.43 


1.53 


1.75 




H{z) (km s^Mpc" 1 ) 


97 


90 


117 


168 


177 


140 


202 




l-cr uncertainty 


±62 


±40 


±23 


±17 


±18 


±14 


±40 





Table 1. Observational H(z) data used to constrain our model. Unless indicated otherwise, the 
values are taken from table 2 in [57]. 

where Oi is the l-cr uncertainty associated with i th point. The observational data used are the 
557 distance modulii of SNela assembled in the Union2 compilation [5] (0.015 < z < 1.40) 
and the 59 distance modulii of GRB from [49] (1.44 < z < 8.10). The combination of these 
two types of data covers a wide range of redshift providing a more complete description of 
the cosmic evolution than the SNela data by themselves. 

3.2 Observational H(z) data (OHD) 

As a function of the redshift, the Hubble parameter is given by 

ff « - -ih% (3 ' 5) 

Therefore it is possible to measure H(z) through a determination of dz/dt. In [55], Jimenez 
et al. demonstrated the feasibility of the method by applying it to a z ~ sample. Simon et 
al. [56], completed by Stern et al. [57] have applied this method using the differential ages 
of passively evolving galaxies to derive a set of 11 observational values for H(z) in the range 
0.1 < z < 1.8. In addition to these values, we also constrain the model using the value of H 
at z = obtained by Riess et al. [58] from the observations of 240 Cepheids, as well those 
obtained by Gazstahaga et al. [59] at z = 0.24, z = 0.34 and z = 0.43 using the radial BAO 
peak scale as a standard ruler (see ref. [59] for more details). The best fit is obtained by 
minimizing the \ 2 function 

2 ro o u l [Hobs(zi) - H th (zi)} 2 
XOHD l^Ao > ^Aoo ' F 0j = 2^ ~2 • ( 3 - 6 ) 



The data used are summarized in table 1. 
3.3 Baryon acoustic oscillation (BAO) 

The use of BAO to test dark energy models is usually made by means of the distance param- 
eter A. However, as pointed out in ref. [47], A is not appropriate to test a model with matter 
production associated with the decays of dark energy. Instead, we can use the dilation scale 



D v (z) 



dz 



H{z) V7 H(z) 



1/3 



(3.7) 



The ratio r s (z c i)/Dv(z), where r s (z^) is the comoving sound horizon size at the drag epoch, 
has been observed at z = 0.35 by SDSS [50] and at z = 0.20 by 2dFGRS [51]. It would be 
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possible to proceed as in refs. [44, 46] and use a fitting formula (developed for the ACDM 
model) to find r s (zd)- However we can avoid this if we follow refs. [43, 45] and minimize the 
X 2 of the ratio Dy Q 35 / Dy 20 given by 

2 r n o i [{ D Va.3s/ D Vo.2a)th ~ (Dv . 35 / Dv . 20) obs] 2 s n 

°0.35/0.20 

In addition to being independent of r s , this ratio is also independent of Hq. The observed 
value for Dy 35 /Dy 020 is 1.736 ±0.065 [51]. 

3.4 Cosmic Microwave Background (CMB) 

The values extracted from the 7- year WMAP data for the acoustic scale (/^(z*)), for the 
CMB shift parameter (R), and for the redshift at the decoupling epoch {z*) can be used 
to constrain the model parameters. The CMB shift parameter and the acoustic scale are 
respectively defined as 

R = \/n mo H 2 r % (3-9) 







and 



^ C z * dz 
71 JO H 
r°° cs dz ' 
J 2* c H 



Ia = • ( 3 - 10 ) 



Since the sound velocity c s is given by 



two additional free parameters are needed to determine the acoustic scale, namely the current 
value of the density parameter of baryons (fi& ) and of radiation (O,^ ). Constraining the 
model with these two additional parameters will require in an increased computational cost. 
However as suggested in ref. [47], we can use the values obtained in the context of the ACDM 
cosmology. This is motivated since the radiation and the baryons are separately conserved, 
and because we want to preserve the spectrum profile as well the nucleosynthesis constraints. 
The observational results from 7-year WMAP data [9] are 

U bo = 2.25 x 10~ 2 h~ 2 and O 7o = 2.469 x 10~ 5 /i~ 2 . (3.12) 

Considering the high redshift values involved here, the dynamical effects of radiation 
cannot be neglected. We therefore must include the density parameter of radiation, Q ro , as 
an extra initial parameter. However, this quantity is related to the density parameter of 
photons through 

7 / 4 \ 3 



n ro = + N cS jn 10 . (3.13) 

Hence, providing that the effective number of neutrino species, N c g, is determined, we still 
have only three parameters to constrain. The departure of N e g from 3 is due to neutrino 
heating by e ± annihilations in early Universe. The value inferred from observations, usually 
close to 3 (3.04 [52], 3.14 [53]), was recently updated to 4.34 [9]. 
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The decoupling epoch occurs when the expansion time becomes less than the Thomson 
scattering time. From a practical point of view, z* may be computed by finding the redshift 
value corresponding to an optical depth of unity 

r(*) = /" TTT^dz = 1, (3.14) 
Jo \X + z)H 

where <jt is the Thompson cross-section. To evaluate this integral, we need to know the 
evolution of the electron density n e as a function of the redshift, which is not simple because 
of its dependence on the recombination process. To avoid to do this complex computation 
each time that one need to know the value of z#, a fitting formula has been developed for the 
ACDM model [54] 

= 1048[l + 0.00124(fi feo /i 2 )-°- 738 ][l+ 5l (O mo /i 2 ) 52 ], (3.15) 

where 

gi = 0.0783^ /i 2 r a238 (l + 39.5(n bo h 2 y°- 763 )-\ (3.16) 
g 2 = 0.560(1 + 21.1(n 6o /i 2 ) 1 - 81 r 1 . (3.17) 

Following refs. [44, 46], we approximate the value of in the A(i)CDM model using eq.( 3.15). 
However, from eq. 3.14, it is obvious that even if we set all the parameters, except Q, to the 
same values, the decoupling redshift for the interacting and the non-interacting cases will 
be different as a consequence of a different redshift-dependence for the Hubble parameter. 
Hence, the validity of this approximation may be questioned; we will come back to this issue 
in the results section. 

Defining v = (I a — R — R obs , z* — z° bs ), the best fit is obtained by minimizing the 
X 2 function 

XcMBpA ,n Aoo ,H ] = vMv\ (3.18) 
where M is the inverse variance-covariance matrix from the 7- year WMAP data 

/ 2.305 29.698 -1.333 \ 
M = 29.698 6825.270 -113.180 (3.19) 
\ -1.333 -113.180 3.414 / 

and the observed value are also taken from the 7- year WMAP data, Ia{ z *) = 302.09 ± 0.76, 
R(z*) = 1.725 ± 0.018 and = 1091.3 ± 0.91. 



4 Results and discussion 

4.1 Case n = 3/2 

The best fit values obtained for the parameters of the ACDM model and the corresponding 
Xmin are shown in table 2 for three different values of N e g : 3.04 , 3.14 , and 4.34. Comparing 
these values to those obtained from the 7-year WMAP data [9] (0 Ao = 0.725 ± 0.016, H = 
70.2 ±1.4 km s^Mpc" 1 ), we find a better agreement when we use N c f[=3.0A. Moreover, the 
Xmin ^ s minimized for this value. For these reasons, we will use iV e g = 3.04 to constrain the 
interacting model. 

For the A(i)CDM model, we present the constraints obtained from each observational 
data set considered separately (distance modulus /i, OHD, BAO, and CMB) as well from two 
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Figure 5. Projection (a) in the plan Oa — , (b) in the plan J7a — -Ho and (c) in the plan f^A^ — Hq 
of the 1-a (darker colour) and the 2-a (lighter colour) confidence region obtained from four different 
types of observational data, the distance modulus fi (blue), BAO (yellow), CMB (pink), OHD (green) 
and the from combination of all of them (C4, yU+OHD+BAO+CMB) (gray). The confidence regions 
obtained from combination of observational constraints C3 (yU+OHD+BAO) are shown in the lower 
panels (dotted contours). The region situated at the right of the dashed line (Qa = 1), involved 
negative matter energy density and must be considered as non-physical. The solid line (f^A^ = 1) 
sets the separation where dark energy is decaying in dark matter (below the line) and where dark 
matter is decaying in dark energy (below the line) . Below the dot-dashed line (Qa = ), k > 
and above, k < 0. The best fit parameters are also indicated by a black dot in each panel. 









3.04 


0.723 
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3.04 


0.718 


c 4 


3.14 


0.712 


C*4 


4.34 


0.635 



fU Aco H (km s _1 Mpc _1 ) 



Xmin 



+0.035 +0.051 1 RO 89 + - 61 + 100 178 

-0.033 -0.051 1 ua.o^ _o,so -0.91 UIU.J.IO 

+0.009 +0.014 -, RQ or +0.49 +0.76 r 7a R oo 

-0.012 -0.019 1 0».00 _o,5g -0.82 OIU.OOO 

+0.009 +0.014 -, RQ 7 r +0.48 +0.75 r 77 n oo 

-0.012 -0.017 1 U». IU _ .57 -0.83 JM.UOO 

+0.012+0.018 1 +0.51+0.78 «1 9 fi q 4 

-0.014 -0.021 1 0».UZ _ 55 _ .83 OlO. 0^1 



Table 2. Best fit values for the ACDM model parameters and the corresponding Xmin f° r three 
different values of N e s inferred from observations 3.04 [52], 3.14 [53], and 4.34 [9]. We have 
considered two different combinations of observational constraint, C3 (/^+OHD+BAO) and C\ 
(/x+OHD+BAO+CMB). The results obtained from C3 arc weakly sensitive to the value of N e g. 
The limits are the 1-a and the 2-a extremal values. 

different combinations of them: C 3 (^+OHD+BAO) and C 4 (/x+OHD+BAO+CMB). The 
best fit values and the corresponding Xmin are shown in in table 3 and the projections of the 
1-a and 2-a confidence regions in the planes £Ia — Qa^ , ^A — ^0 an d ^A^ — Hq , are shown 
in figure 5. 

We consider the combination of constraints C3, which excludes the CMB constraints, 
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Figure 6. Percentage of difference on the value of the integral J** dz/(H(z)(l + z)) evaluated for 
the ACDM model and for the A(t)CDM model for a corresponding value of 51a (z* is given by 
eq. (3.15)). A positive difference means that the ACDM value is larger than the A(t)CDM one and a 
negative difference, the converse. The dashed lines and the solid lines represent respectively, the 1-a 
and the 2-a confidence contours obtained from the CMB data, and from the combination of data C4 
(/i+OHD+BAO+CMB). 



Oa 



Ha 



Ho (km s^Mper 1 ) 



Xmin 



xi 



xi 



0.746 
0.804 
0.812 
0.710 
0.753 
0.720 



+0.178 
-0.122 

+0.196 
-0.570 

+0.188 
-0.775 

+0.115 
-0.189 

+0.142 
-0.106 

+0.013 
-0.013 



+0.254 
-0.172 

+0.196 
-0.690 

+0.188 
-0.812 

+0.161 
-0.315 

+0.229 
-0.153 

+0.020 
-0.021 



+0.770 +01.383 
-0.355 -00.461 



+0.970 +01.6CO 
-0.970 -01.090 



1.061 
1.174 

0. 681 

1.011 

1.099 
1 nnfi +0.020 +00.030 

1. UUD _o,019 -00.029 



+2.743 +40.695 
-0.675 -00.680 

+0.254 +01.116 
-0.038 -00.046 

+0.549 +00.975 
-0.301 -00.408 



69.94 
72.68 

68.89 
69.97 
69.77 



+00.97 +01.49 

-00.87 -01.30 



+07.99 +10.50 
-11.67 -16.06 



+ 14.08 +23.80 
-11.42 -16.11 

+00.91 +01.38 
-00.84 -01.25 

+00.52 +00.79 
-00.52 -00.79 



565.834 
7.970 
0.000 
0.015 
575.888 
576.267 



565.878 566.182 
8.824 8.790 
1.185 1.240 
0.054 



Table 3. Best fit values for the interacting model parameters and the corresponding Xmin 
associated with each data set and two combinations of them: C3 (//+OHD+BAO) and C4 
(/i+OHD+BAO+CMB). The two last columns, xi an d xli represent the partial contribution of 
each data set to the value of C3 and C4, respectively. For the sake of comparison, the values obtained 
for the ACDM model using the same constraints are also shown. The limits are the 1-a and the 2-a 
extremal values excluding the non-physical region J7a > 1- 

because there is a possible circularity problem in our analysis of the CMB data. Indeed, we 
have used the values of £lb Q , ^70 an d Q ro obtained in the context of the ACDM model and we 
have approximated the value of using a fitting formula developed for the same model. In 
the latter case, the difficulty to obtain the exact value of comes from the presence of the 
electron density n e , whose value at different redshifts depends on the recombination history, 
in the computation of the optical depth (eq. (3.14)). Hence, in order to have a rough idea 
of the error induced by this approximation, we have simply removed the electron density 
from eq. (3.14) to get the integral J^* dz/(H(z)(l + z)) (where z* is given by eq. (3.15)) 
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and compared the values obtained for the interacting and the non-interacting cases. If the 
difference is not too large, that means that up to z#, the evolution of the Hubble term and thus 
of the scale factor, are similar in both cases; this in its turn implies that the recombination 
history has also to be similar since for the comparison, all the initial parameters are fixed to 
the same value (except Q). In figure 6, we see that for the CMB constraints considered alone, 
the percentage of difference reaches up to ~ 7.5% in the 2-cr confidence region and ~ 5% in 
the l-o" confidence region. When all the data are considered, the percentage of difference in 
the 2-cr and in the 1-a confidence region falls respectively below 2.5% and 2%. Hence, we can 
expect that eq. (3.15) provides a reasonable approximation for z#, at least in the confidence 
regions obtained from the combination of constraints C4. 

We can notice that the confidence regions obtained from BAO (l-o") , from OHD (l-o - ) 
and from the distance modulus p (2-a) extend partially beyond = 1- However, a value 
of f^Ao greater than 1 would imply a negative value for the current energy density of matter 
(p mo < 0) and therefore must be considered as non-physical. Similarly, we can also notice 
that the l-o" confidence regions obtained from BAO and OHD extend (marginally in the case 
of OHD) into region of the parameter space where k > (below the dot-dashed line on figure 
5). We have shown in section 2.1 that the solutions with k > necessarily reach a point 
where p m becomes negative if we are looking sufficiently far in the past. However, for the 
redshifts involved in the computation of the BAO and the OHD constraints, this point is 
never reached and p m remains always positive. Hence, in this case, there is no reason to 
consider these constraints as non-physical. In any case, these two regions are excluded by at 
least 2-o" from the combinations of constraints C3 and C4. 

Regarding the best fit values obtained, the processes of decay of dark matter into dark 
energy (^Aoo > 1) is favoured by only one constraint (BAO), while the decay of dark energy 
into dark matter (^a^ < 1) is favoured by the three other, as well as by the combinations 
C3 and C4. However, all the l-o" confidence regions extend below and above the the plane 
^Aoo = 1 (which corresponds to the ACDM model) and therefore both processes still are 
statistically allowed. 

As shown in section 2.1, we need to have S]^ ~ £Ia in order to relate and explain 
the cosmological constant problem and the coincidence problem. However, as we can see on 
figure 5, these values are excluded by at least 2-ct due to constraints C4. Actually, the closest 
point to the line SIa^ = ^A situated on the 2-0" contour is Oa = 0.718, ^Aoo = 0.983. 
For these parameters, the value of pa remains mostly constant (paJ PA q ~ 1-04), and a 
coincidence problem is still present since the current ratio of matter and dark energy density 
(,Pm /pAa = 0.393) lies in the transition zone between the point where dark energy starts to 
dominate and the late-time phase where the ratio p m /pA becomes approximately constant 

{Pm^/PK^ =0.017). 

The best fit values obtained for the ACDM model from the constraints C4 (table 2) are 
included in the l-o" confidence region of the interacting model. Moreover, the Xmin ^® only 
slightly smaller for interacting model (576.267 vs. 576.638). Hence, we may wonder if the 
introduction of an extra-parameter (^a^) is really justified. The improvement of the Xmin 
value may be assessed by the mean of the Bayesian information criterion [60] , defined as 

BIC = -2 In £ max + K In N = ( X 2 min + C) + K In N, (4.1) 

where C max is the maximum likelihood, K is the number of parameters for the model (2, 
for the ACDM model, 3 for the interacting model), N the number of data points used in 
the fit (N = 635) and C a constant independent of the model used. Following ref. [61], 
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we will regard a difference of 2 for the BIC as a non significant, and of 6 or more as very 
non-significant improvement of the x min value. Since we get Abic = 6.08, we can conclude 
that the addition of an extra parameter is not warranted by the marginal decrease in the 
value of Xmin- 

If we consider now the combination of observational constraints C3 (which excludes the 
CMB data), the improvement of the x m in va hie remains insignificant (Abic=6-15) and the 
closest point to the line = &a on the 2-a contour is now Oa = 0.615 , Qa^ = 0.693, 
which leads to the ratio pa J p\ Q ~ 3.37, still far to being able to explain the discrepancy with 
the theoretical value. The current value of the ratio p m /pA is now closer to the late-time 
value (0.443 vs. 0.626), but not sufficiently to explain the coincidence problem. 

4.2 General case 

The motivation to study the model with n = 3/2 was purely phenomenological and did not 
rely upon any deep principles. Then, it will be interesting to extend our analysis to general 
values of n. However, we do not perform an analysis as detailed as for n = 3/2 and restrict 
it to constrain the model parameters in the range n = to n = 10. The results are presented 
in figure 7 for the combination of observational constraints C3 and C4. 

We can see on figure 7 that the best fit values for Qa an d Hq are nearly constant for 
both C3 and C4, and for the latter, this is even the case for the confidence intervals. These 
results seem a priori very surprising, but could be understand if the constant Q is sufficiently 
small (though still non-zero). Indeed, in the limit of a small Q, it is possible to reproduce 
exactly the same cosmological evolution, for t < to, using two different values of n. 

If we suppose that the term Qp\ can be neglected in the matter energy conservation 
equation (eq. (2.2)), we find a dependence on n only in the energy density p\ and in the 
quantities derived from it. We are then interested to determined the conditions needed to 
have PA n {t) = PA A (t), with n ^ n. In the limit where Q is small, the solutions found for pa 
(eq. (2.3)) become 

pAn(t) = PAo-Qpl t- ( 4 - 2 ) 
Hence the condition to have pA n (t) = pa a (t) is given by 

Q = QP17' ( 4 - 3 ) 

which provides an explanation to the exponential relationship observed on the lower panels 
of figure 7 between the best fit value of Q and n. 

For the best-fit values, the same cosmological evolution is closely reproduced indepen- 
dently of the value of n, hence we find the same value for x m hr Concerning the 1-a and 
2-a confidence regions, the results are weakly sensitive to n for the constraints C4, (for the 
higher values of n, Q becomes sufficiently large to start to observe a weak deviation from 
the exponential relationship) but are more much important for the constraints C3. However, 
even in this case, the constraints obtained do not allow pa to decay from an initially large 
value to that observed today. 

5 Conclusion 

In this paper, motivated by the hope to shed some light on the coincidence and the cos- 
mological constant problems, we have applied the phenomenological decay law p\ = —Qp A , 
originally proposed in ref. [19] to described the decay of dark energy into radiation, to study 
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(a) (b) 




Figure 7. The panels in the first three rows show the best fit values (solid line), the \-a (dark shade) 
and the 2-cr (light shade) confidence intervals for fiA , Hq and Q/Qbfv at different value of n. In 
the panels of the last row, the dots represents the best fit values Q b fv & n d the solid line is given by 
eq. (4.3) where we have used the best fit value found in section 4.1 (h = 3/2) for Q and p\ . For 
Q/Qbfv, since Qbfv turns out to be always negative, Q/Qbfv > corresponds to dark matter 
decay and Q/Qbfv < to dark energy decay. These results were obtained using the combination of 
constraints C 3 (/i+OHD+BAO) in (a), and the combination of constraints C 4 (^+OHD+BAO+CMB) 
in (b). 

the interaction between the dark sectors. In order to ameliorate these two problems, we were 
primarily interested in the decay of dark energy into dark matter (Q > 0) , but we have also 
considered the decay of dark matter into dark energy (Q < 0). 

From dimensional analysis, we have shown that the model with n = 3/2 (and only 
this model) leads to solution where the ratio of energy densities p m /pA is constant at late 
time. An important feature of this model is the possibility to have, in addition to this late 
phase (which could explain the coincidence problem), an early phase qualitatively similar to 
what is obtained in the ACDM model (which already provides a good fit to the observational 
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constraints). 

To constrain the three free parameters of the model with n = 3/2, i.e. the current 
density parameter of dark energy (O,\ ), the late-time density parameter of dark energy 
(QAoo — related to Q) and the current Hubble parameter (Ho), we have mainly followed 
the procedure of refs. [43-46] and considered the observational constraints involving the 
distance modulus p of type la supernova (SNela) and gamma-ray bursts (GRB), the baryon 
acoustic oscillation (BAO), the cosmic microwave background (CMB) and the observational 
Hubble rate (OHD). The constraints obtained from the CMB data could be affected by a 
circularity problem. To constrain the A(i)CDM model with the CMB data, we have used 
some parameters obtained from the usual ACDM model (^6 , 70 , £l ro , Although we 
have argued that it was legitimate to do so, that could possibly lead to biased results. For this 
reason, we have presented our result considering two different combinations of observational 
constraints, C 3 (^+OHD+BAO) and C 4 (/i+OHD+BAO+CMB). 

For the constraints C4, the best fit parameters are given by Qa = 0.720, SIa^ = 1-006 
and Hq = 69.77 km s _1 Mpc _1 . Since f^A^ > 1 corresponds to Q < and f^A^ < 1 to 
Q > 0, the process involved at the best fit point is the decay of dark matter into dark energy 
and in fact exacerbates the coincidence and the cosmological constant problems. We have 
shown that in order to explain both the coincidence and the cosmological constant problems, 
we need to have £l\ ~ ^Aoo- The point lying in the 2-a confidence region where these 
parameters are closest is Oa = 0.718, Oa^ = 0.983. For these values, p\ remains nearly 
constant [phj Ph ~ 1.04) and the current value of the ratio p m /pA is far from the late-time 
value (0.393 vs. 0.017). 

For the constraints C3, the best fit parameters (^a = 0.753, ^Aoo = 1-099 and Hq = 
69.97 km s _1 Mpc _1 ) also correspond to the decay of dark matter into dark energy. The 
closest point (in the 2-a confidence region) to the line Q\ = SIa^ is now Q\ = 0.615, 
OAoo = 0.693. At this point, the value of PA;/PA remains far away from the desired value 
(sa 3.37), but the current and late-time value of the ratio p m /pA are now of the same order 
of magnitude (0.626 and 0.443, respectively). However, to provide a convincing explanation 
to the coincidence problem, we should have p mo /pA ~ Pm^l Woo) which is not the case here. 

From a statistical point of view, we do not find that the addition of an extra parameter 
was justified. The x m in va l ues obtained for the A(i)CDM model were only slightly smaller 
than the value obtained for the ACDM, 575.888 vs. 576.178 for C 3 and 576.267 vs. 576.638 for 
C4. In both cases, according to the Bayesian information criterion (BIC), the improvement 
of the Xmin i s considered to be insignificant (Abic = 6.15 for C3 and Abic = 6.08 for C4). 

Finally, we have extended our results to generic values of n ranging form to 10. The 
values of x m in an( ^ the corresponding values for Vl\ and for Hq are nearly constant, whereas 
the value of Q have a exponential dependence on n. We have shown that these results 
may be understood in the limit of a small Q. Since the same values of x m in are f° un d, the 
conclusion that we drew for the BIC holds for any value of n, i.e. the improvement of the 
Xmin i s insignificant. Moreover, even if we consider the 2<7-confidence region, the constraints 
obtained do not offer any explanation for the cosmological constant problem. 

Since when we consider the exclusion limits obtained from observational data, the model 
studied in this article fails to provide an explanation to the coincidence and the cosmological 
constant problems, and since this model is also disfavoured by the BIC, we can conclude that 
the usual ACDM model remains the most reasonable description of the Universe. 
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